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SUMMARY 


Turbulent flow in a duct in which the cross-section undergoes a transi- 
tion from square to round has been measured at a Reynolds number of 35 350. 

The hydraulic diameter, 40 mm, was the same in the square inlet and round exit 
sections, resulting in a cross-sectional area reduction of 21.5^. The tran- 
sition takes place over two hydraulic diameters. Laser Doppler velocimetry 
was used to measure the streamwise and secondary mean velocity components 
and the associated turbulence levels and cross-correlations in stations up- 
stream of, in, and downstream of the transition. The results provide an un- 
derstanding of the flow and assist the evaluation of numerical techniques. 

The boundary layers at inlet and exit of the transition were approxima- 
tely 13% and 20 % respectively of the hydraulic diameter. Thickening of the 
boundary layers at the duct corners is associated with the secondary flow 
away from the corner fillets and along the duct walls, towards the duct mid- 
planes, Cross -stream velocities of up to 7 % of the bulk velocity were mea- 
sured: the secondary flow is produced by the longitudinal curvature of the 
wall and the related pressure gradients. 



1. INTRODUCTION 


Flows in ducts which change cross-sectional area and shape are used in 
a wide range of engineering configurations and this report is concerned with 
the flow in a transition from a square to a circular cross-section with a de- 
crease in area of 21.5% over two hydraulic diameters. The lateral convergence 
of the flow is accompanied by the generation of cross-stream velocities, or 
secondary motions, which result from the pressure fields 

In long ducts of square section, secondary flows of the "second kind" 
arise and are driven by gradients of Reynolds-stresses with a magnitude of 
less than 1.5% of the bulk flow velocity: see, for example, reference 1. 

In the present flow the turbulent boundary layers at the inlet to the tran- 
sition piece are thin and it can be expected that any secondary flows of the 
second kind will be much less than, say, 1% of the bulk flow velocity. Any 
larger cross-stream velocity within the transition piece stems from the 
change in cross-sectional shape and the related cross-stream pressure gradients. 

The measurements of reference 3 provide some information of the flow 
in a transition from a rectangular to a circular section (and vice versa) 
with no change in cross-sectional area. The data includes contours of 
streamwise isotachs and the streamwise pressure variation. Two configura- 
tions were investigated in which the transition occurred over 0.69 and 
2.76 hydraulic diameters and it is evident that the length of the transi- 
tion is influential on the flow development. For example, cross-stream ve- 
locities presented for the round-to-rectangular transitions reached maxima 
of about 0.18 and 0.09 of the maximum streamwise velocity for the short and 
long transition ducts respectively. For the rectangular entry sections, 
the flows were fully developed at the inlet to the transition ducts and 
thus stress-driven secondary flows were also present. 

The flow investigated in this report is relevant to the case of air- 
craft intake ducts. The measurements are intended to quantify the mean and 
fluctuating components of the streamwise and cross -stream velocities of the 
flow and so provide better understanding and a basis for the testing of de- 
sign methods. In common with the previous investigations of references 1 
and 2, laser Doppler velocimetry was used to measure the velocity components. 

The following section describes the experimental procedure and section 
3 presents the measurements. Concluding remarks are provided in section 4. 


2. FLOW CONFIGURATION AND INSTRUMENTATION 

The transition from a 40 mm square cross-section to a 40 mm diameter 
circular cross-section took place in a duct length of 80 mm. The duct is 
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depicted in figure 1 together with the coordinate system adopted in the re- 
port. Each cross-section in the transition is formed by the intersection 
of a square and a circle, which resembles, to within ^ .S% of the radius 
a super ellipse with a shape factor of unity. The radius of the circle de- 
creases linearly along the transition from 28.28 mm to 20.00 mm. The tran- 
sition piece was installed in the water tunnel previously used in reference 2. 
The bulk velocity through the square section was 0.885 m/s which corresponds 
to a Reynolds number of 35 350, based on the hydraulic diameter. The tempe- 
rature of the water was maintained at 20 t 2° C for all the experiments. 
Measurements of velocity were restricted by limits of optical access, caused 
by refraction at curved surfaces within the transition piece and in the 
downstream pipe. 

The optical arrangement of the laser Doppler velocimeter was identical 
to that of reference 2. Doppler signals were processed by a frequency- 
tracking demodulator (Cambridge Consultants CCOl). A careful appraisal of 
the measurement errors has already been made in reference 2, where the preci- 
sion (random error) of the mean and fluctuating components of velocity was 
assessed as t 1J% and from t 1J% to ± 5% respectively. In interpreting the 
results of the next section, it is important to note that the accuracy (sys- 
tematic error) of the secondary velocity components is limited by the align- 
ment of the velocimeter to the duct (see reference 1). The resulting syste- 
matic error is not more than 0.008 V . 


3. RESULTS A>iD DISCUSSION 

Symmetry tests on the streamwise flow in the duct upstream of the tran- 
sition piece were carried out, with the result that measurements at equiva- 
lent locations on either side of the symmetry plane of the duct agreed with- 
in the precision of the measurement. 

Figure 1 shows the development of the streamwise velocity on the centre- 
line of the duct. The velocity increases from 1.12 upstream of the tran- 
sition to 1.43 downstream, as a consequence of the decrease in cross-section. 
Figures 2a-2e show profiles of the streamwise mean velocity, U, and the cor- 
responding root mean square velocity, u, at a number of stations upstream, 
within and downstream of the transition piece. The boundary layer thickness 
on the symmetry plane, defined at 95% of the maximum velocity, is about 13% 
of the hydraulic diameter at the inlet to the transition piece (figure 2a), 
but rises to about 20% at the exit (figure 2d). The change in the profiles 
is comparatively small, as would be expected from the centreline change in 
figure 1. Careful examination of the profiles near the sidewalls shows 
that the fluid near the planes xp/D and X3/D = 0 accelerates relative to that 
near to the fillets which form the transition, and that the boundary layers 
are thinner away from the fillets. The profiles of Q (and v, figure 2f) are 
those expected for a turbulent wall boundary layer. 
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Figure 3 presents profiles of the cross -stream velocity for three sta- 
tions within the transition piece. The measurements quantify the symmetry 
about planes at 45® to the X£ or x, planes, which is within the limits of 
the systematic error discussed in'^section 2, as are the magnitudes of the 
cross-stream velocities on the symnetry planes (x 2 /D or X3/D=0). The mea- 
surements show that large values of the cross-stream velocity are found 
close to the walls and to the fillets which form the transition from the 
square to the circular cross-section. The largest values shown are of the 
order of 0.07 V^., although the shape of the profile suggests that the maxi- 
mum lies closer to the wall than X3/D=0.9. The magnitude and direction of 
these velocities indicates that their source is pressure-gradient - rather 
than Reynolds stress-gradient. The direction of the cross-stream velocities 
is always away from the channel diagonal (line A-A' on figure 4) and thus 
cannot satisfy continuity in the cross-stream plane. This implies an acce- 
leration of the streamwise velocity component near the symmetry planes (li- 
nes B-B' on figure 4) which is indeed observed, particularly in figure 2(c). 

The origin of the lateral pressure-gradient, which drives the cross- 
stream flow, is due to the increase in pressure near the fillet relative to 
that in the region of the symmetry planes. This increase is brought about 
by the streamline curvature occurring in planes containing the channel dia- 
gonal and the duct centre-line upon entry to the transition piece, as shown 
in fig. 4 (a). In contrast, no such curvature occurs with streamlines con- 
tained in planes such as that depicted in Fig. 4 (b). It is noted that this 
mechanism should be reversed on exit from the transition piece and can be 
expected to cause the rapid decay of the cross-stream components. 

These lateral pressure variations were also observed by Mayer (3); the 
qualitative changes in pressure and the related secondary flow directions 
are in accordance with the above discussion. However the simultaneous oc- 
currence of longitudinal convergence on one wall and divergence on the other 
wall of the transition duct results in cross-stream flows not directly com- 
parable to those presented in this report. 

An order-of-magnitude estimation of the affect of curvature, given in 
Appendix 2, suggests secondary flow velocities not more than 0.16 of the 
bulk velocity, is, considering the crude nature of the analysis, in 

good agreement with experiment. 

Measurements of the uw and uv cross correlations were also obtained. 

As expected they were very low (< 0.2 x 10"3v^.2) over the central region of 
the flow, rising to about 0.5 x 10"3Vp2 pear the sidewalls. These values 
are close to the limits of accuracy of the velocimeter and are thus not pre- 
sented. 

The results from which the figures have been plotted are tabulated in 
Appendix 1 . 
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4. CONCLUDING REMARKS 


1. The experiments have quantified the streamwise and cross-stream ve- 
locities in a square to circular transition piece with thin inlet boundary 
layers. 

2. The streamwise boundary layers assume smaller thicknesses away from 
the fillets which form the transition piece. The cross-stream velocity 
profiles are consistent with this observation; magnitudes of = 0.07 of the 
bulk velocity have been measured. 

3. The secondary velocities are driven by cross-stream pressure gradients 
induced by the longitudinal curvature of the streamlines on passing through 
the transition piece. 

4. Because of the importance of wall curvature, its accurate representa- 
tion is crucial to the calculation of the flow. 


5 



APPENDIX 1 


TABULATED DATA 


TABLE I - SQUARE-TO-ROUND TRANSITION MEAN VELOCITY, 
TURBULENCE LEVEL MEASUREMENTS 

REYNOLDS NUMBER 35,350 

(All quantities normalised by the bulk 
velocity, V^ = 0.885 °^/5) 



D 

m 

Hi 

■ 

u/V 

c 

xlO^ 

■ 

a/v^ 

u/Vc 

- 

tr/v^ 

— 

u/V^ 

xlO^ 

-2.0 

0.0 

0.0 

1.118 

0.721 

-1.0 

1.136 

0.721 

0.0 

1.145 

0.722 



0.1 

1.118 

0.721 


1.136 

0.721 


1.145 

0.768 



0.2 

1.118 

0.721 


1.136 

0.721 


1.145 

0.811 



0.3 

1.118 

0.721 


1.136 

0.721 


1.145 

0.839 



0.4 

1.118 

0.721 


1.136 

0.811 


1.145 

0.901 



0,5 

1.118 

0.946 


1.133 

1.117 


1.138 

1.533 



0.6 

1,107 

2.389 


1.116 

2.636 


1.109 

3 .426 



0.7 

1.050 

4.704 


1.068 

5.273 


1.025 

5 .409 



0.8 

0.964 

6.200 


0.964 

6.698 


0.960 

6.220 



0.9 

0.865 

7.341 


0.879 

7.910 


0,865 

7.697 

3.0 

0,0 

n 

1.421 

0 .914 

4.0 

1.430 

0,959 







1.421 

0.914 


1.430 

0-959 







1.421 

0.959 


1 432 

1.005 






IB 

1.421 

1.074 


1.435 

1.188 







1.414 

1.736 


1.427 

1 .873 






0.5 

1.391 

3.290 


1.400 

3 .610 






0.6 

1.343 

4.769 


1.352 

4.769 






0.7 

1.288 

5.636 


1.297 

5.636 






0.8 

1.220 

6.214 


1.220 

6.503 






0.9 

1.119 

7.949 


1.087 

9.683 
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0.0 
- 0.1 
-0 . 2 
-0.3 
-0.4 
-0.5 
- 0.6 
-0.7 
- 0.8 
-0.9 
-0,95 


1.172 

1.172 

1.174 

1.179 

1.172 

1.167 


0.960 

0.842 


1.255 

1.260 


1.246 

1.217 

1.145 


X 

x1q2 

o < 
too 

xlo^ 

0.766 

-0.646 

0.572 

0.811 

.646 

0.682 

0.766 

-0.646 

0.706 

0.856 

-0.807 

0.725 

0.991 

-0.646 

0.862 

1.690 

-0.646 

1.646 

3.516 

-0.807 

2.433 

5.341 

-0.646 

4.327 

6.693 

-0.807 

5.327 

6.761 

-0.807 

- 

0.901 

-0.646 

0.745 

0.901 

-0.606 

0.719 

0.946 

-0.646 

0.765 

1.127 

-0.646 

0.778 

1.149 

-0.646 

1.098 

2.118 

-0.646 

1.698 

3.516 

-1.055 

3.568 

5.612 

-0.969 

4.301 

6.761 

-2.642 

5.409 

7.843 

-3.232 

5.557 

1.680 

-O . 646 

2 .042 

1.805 

-0.646 

1.805 

2.251 

-0.646 

1.953 

2.299 

-0.646 

2-150 

2-411 

-0.969 

1.850 

3.074 

-1.131 

2.444 

4.395 

-1-707 

3.280 

5.950 

-2.262 

4.316 

7.234 

-3.706 

5.334 

7.978 

-5.171 

5.601 

0.811 

-0.323 

0.626 

0.811 

-0 . 646 

0.626 

0.811 

-0. 646 

0.626 

0.901 

-0.323 

0.609 

1.262 

-0.646 

0.691 

1.848 

-0 . 646 

2.307 

3.651 

-0.484 

3 .700 

5.544 

-0.807 

4.533 

6.558 

-0.807 

5.648 

7.572 

-0.807 

6.154 

11.630 



-0.807 

5.475 

1 — 



-0.637 

-0.637 

-0.637 

-0.637 

-0.797 

-0.797 

-0.637 

-0.797 

-0.637 

-0.637 


-0.637 

-0.637 

-0.637 

-0.637 

-0.637 

-0.637 

0.000 

0.000 

1.530 

3.188 


-0.159 
-0.159 
0.000 
-0.159 
-0.159 
0.556 
0.797 
2.231 
3 .188 
5.206 


-0.318 

-0.318 

-0.318 

-0.318 

-0.318 

-0.318 

-0.797 

-0.797 

-0.797 

-0.797 

-0.797 


0.867 

0.867 

0.867 

0.867 

1.068 

2.134 

3.584 

4.886 

5.288 

5,685 


0.819 

0.819 

0.959 

0.959 

1.066 

1.753 

3.231 

5.048 

5.171 

5.546 


2.050 

2.203 

2.209 

2.475 

2.147 

2.721 

3.186 

4.472 

5.613 

6.168 


0.819 
1.005 
0.968 
1.096 
1.230 
2.303 
3 .778 
5.026 
5.719 
2,689 


*For the w and w measiirements the absolute values of the coordinate x 
correspond to -x^ and of the coordinate /_ to x^/„ . 
































^2/d 



u/v 

u/V^ 

xl02 

V/V^ 

xlO^ 




1.258 

1.258 

1.258 

0.856 

0.856 

0,901 

-0.646 

-0.646 

-0.646 


1.189 

1.113 

1.036 


1.307 

1.307 

1.307 

1.307 

1.307 

1.298 

1.256 

1.197 

1.121 

1.030 


v/v^ 

xlO^ 

Wv^ 

1 xlO^ 

0.672 

1-0.318 

0.699 

‘-0.318 

0.719 

-0.159 

0.765 

-0.159 

! 1.142 
1 

-0.159 



0.901 

0.901 

0.946 

0.991 

1.442 

3.110 

4.703 

5.986 

6.699 

7.269 


0.856 

0.901 

0.901 

0.991 

1.127 

2.402 

4.192 

5.544 

7.572 

8.056 


1.051 

1.051 

1.033 

1.051 

1.678 

2.369 

3.698 

5.175 

5.367 

5.803 



-0.646 

0.719 

-0.478 

0.872 

-0.646 

0.745 

-0.478 

0.961 

-0.646 

0.773 

-0.318 

0.959 

— 0 . 646 

0.786 

-0.318 

1.005 

-0.646 

1.210 

-0.159 

1.883 

-0.807 

1.743 

0.000 

2.248 

-1.417 

3.542 

0.318 

4.235 

-2.101 

4.955 

1.274 

5.482 

-3.210 

5.487 

2.039 

5.589 

-3.393 

5.784 

3.505 

6.139 
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APPENDIX 2 


Estimation of the cross-stream pressure variation 


Referring to figure 5, consider the centripetal force balance due to the stream- 
line curvature for the fluid particle shown : 

6 m ^ = <5p ( 6 x*<Sy) (1 ) 


{Nhere 


Hence 


6 m = p( 6 x* 6 y 6 r) is the mass of the fluid particle, 

6 p is the difference in pressure across the particle in the r direction 
R is the radius of curvature of the streamline} 


^ 6 r 

Using the following approximations 


( 2 ) 


V « U 


B 


{where 


Ug is a representative streamwise velocity, 
D is the hydraulic diameter of the duct 
P 2 and p-| are defined in figure 5} 

D 

Pe'Pi “ " it I 


(3) 


Considering motion in the cross -stream plane (refer to figure 4 (c)) and 
applying Bernoulli's equation between points near the fillet and points near 


X 2 = X 3 - D 


aU - 






{where aU is the induced cross -stream velocity} 

Eliminating the term for (P 2 “P-|) in (4) using (3) there results 

/ II- II 

AU 


AU 

'Ji 




for D = 0.04 m. 


An estimate for R, derived in figure 5, is R > 1*5 m. and hence ^ < 0.16: 
Viscous forces will act so as to lessen the magnitude further 
still. 


aU 

Ob 


(4) 


(5) 
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APPENDIX 3 - NOMENCLATURE 


defined as Dj^/2 
hydraulic diameter (40mm) 

mean component along X-j co-ordinate direction 
r.m.s. component along X-j co-ordinate direction 
cross correlation of fluctuation in X-|'''X 2 plane 
cross correlation of fluctuation in X^'^^X 2 plane 
mean component along X 2 co-ordinate direction 
bulk mean velocity defined in square section 
r.m.s. component along X 2 co-ordinate direction 
mean component along X^ co-ordinate direction 
r.m.s. component along X^ co-ordinate direction 


defined in figure 1 
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x^/D 



(a^ Cb) Cih* I.S (o') D^s 1.75 


Fig. 3 Secondary velocity profiles 



(a) Section AA' 


3 ^ 


C 


(b) Section BB' 

(not to scale) 

Fig. 4 Effect of longitudinal streamline curvature on cross-stream pressure variation 


r 



Note: radius of curvature of streamline approximated as radius of circle which 
inscribes points A,B,C. Hence Rsl.5 m. 

Fig. 5 Definition sketch for the estimation of the cross-stream pressure variation 
(Appendix 2). 
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